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INTRODUCTION 



This study describes the development of a hybrid program for 
obtaining **3-dimensional'* surfaces of temperature as a function 
of time and distance in an air-cooled glass slab. The purpose of 
the program is to demonstrate an effective use of the TR-48/DES-30 
Desk- Top Analog/Hybrid Computing System as applied to the simu- 
lation of a system described by a partial differential equation and 
towards the improvement of the man-machine interface by providing 
a visual, "3 -dimensional" display. 

The logic for arranging the solutions from the TR-48 General 
Purpose Analog Computer in **3-D" format, and for controlling the 
plotting process automatically is supplied by the DES-30 Digital 
Expansion System. This general -purpose, low-cost digital logic 
package, especially designed for operation with desk-top analog 
computers, also provides mode control and other logic tasks as 
required in the high speed repetitive operation (rep op) solution of 
the system equations. 

The DES-30 System, although capable of operating autonomously for 
use in digital design or instruction, was designed primarily as an 
expansion to the TR-48 analog computer to provide basic hybrid 
capabilities to the small computer facility. It also can be combined 
readily with other general-purpose digital computers to provide 
extended control functions. Thus, it makes available to the analog 
computer in this application, the necessary digital logic functions to 
implement this essentially hybrid program. 
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SYSTEM DESCRIPTION 



and equation (1) reduces to 



The physical system of interest is shown in Figure 
1. The problem is to investigate the cooling (by air) 
of a slab of glass (subsequent to heat treatment) . 
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Figure 1: Glass Slab Being Air Cooled 

In particular, it is necessary to determine how cold 
the ambient air can be (how low T^ can be) before 
severe temperature gradients in both space (x) and 
time (t) are developed. (Extreme gradients result 
in stress cracks). For this, it will be desirable 
to observe a surface'* oftemperatureTasa func- 
tion of distance x and time t. 

Cooling is described by the one-dimensional heat 
diffusion equation: 



01 3t 



(1) 



where 



T = f(t,x) 
k 



a = 



pc 



= specific heat 

= density 

= conductivity 



METHOD OF SOLUTION 



Normalizing the Variables: Equation (1) can be Sim- 
plified by using dimensionless variables as follows: 



3T 

be 



(4) 



Obfaining T{0) and T(z}: PDE(4) is Solved in two dif- 
ferent ways simultaneously: first, it is solved by 
finite differencing the z dimension (corresponding 
to distance); there results a set of basic ally identical 
equations (see Appendix I for derivation of the equa- 
tions at the air-glass interface): 

dT. T - 2T. +T. 

= -ii J lli(j = 0,1,2,3) (5) 

(Az)' 

where j is the number of a given point in z. 
Refer to Figure 2. 
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Figure 2: Definition of Points for Space Finite Differencing 

Second, it is solved by finite differencing the 6 
dimension (corresponding to time) to get the set of 
equations: 



d^T, 



dz 



T - T 
_i izl 

Ae 



(i = 1,2,3,4)* 



(6) 



where i is a given instant in time. 



Z = — (normalized distance) 

1j 

(2) 

9 = t (normalized time) 

This results in 

2 2 

a T 1 a T 

2 2 2 
Sx L Bz 

(3) 

BT _ g dT 

dt ^2 ae 



The solution of equations (5) (''parallel" solution) 
are T j {$) which describe temperatures as a fimction 
of time for a set of discrete points in distance (the 
index j identifies the point). 

Equations (6) (''pseudo serial" solution) yield 
Tj^(z) which describe temperature as a function 
of distance for a set of discrete points in time 
(the index i identifies the point). Note that, for 
this solution, computer time plays the role of dis- 
tance , z. While all T^iz) are produced simultan- 
eously with respect to computer time, they really 



*/ sfarts wifh .7 rather tfion 0 since at time 6= 0 (t = 0) the slab is at 
° known initial temperature, (z). 
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represent snapshots" of temperature distribution 
in the slab at progressively later points in real 
time, 

A simplified implementation of one equation each 
of sets (5) and (6) is shown in Figure 3. 
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Figure 3: Typical (Simplified) Analog Blocks 

Note from the computer diagram that the serial 
solution is basically unstable. If it runs too long 
(exceed the distance z =1), it will surely blow 
up (diverge). It will also tend to magnify noise. 
This is one of the major reasons the pseudo- serial 
solution is rarely used. 

Inifial Conditions: Another reason for not using the 
pseudo-serial solution, is that the I.C.'s required 
for it are usually difficult to come by. In this case, 
the LC.'s are: 



T. (z = 0) 



for all i (all time) 



dz 



(z = 0) for all i (all time) 



Because of the symmetry of the system, it can be 
said that: 



dT. 



(since, at the center of the slab-r-which is the 
center of symmetry — the rate of change of 
temperature with respect to distance is zero). But 
Ti(z =0) is the temperature atz = 0 (center of slab) 

at time to, t^, tj, 1^ (where n = 4, in this case). 

We can obtain these LC/s from the parallel solution by having 
track/stores pick up the values T^(d) at the instant ti and de- 
liver these values to the LC, inputs of the second integrator 
[Tj^iz)] in each block of the serial solution. See Figure 4 
below. 




^ (- = o)=o 



Figure 4: Program to Obtain I.C. for the Semi Serial Solution 

The DES-30 System generates the proper TRACK- 
HOLD signals for these track/ stores during each 
operate (OP) period. At four equally spaced in- 
stances during the OP period, a monostable timer 
(MT) is set for about 100 /xs (the time necessary to 
assure proper tracking). Each time this happens, a 
4-bit ring shift register advances one position. 
The position of the circulating bit in the register 
is used to sequence the MT signal to the four track/ 
stores. 

The I.C.'s for the parallel solution are Tj (0=0), 
the temperature distribution in the slab at time 
0=0, which is known, 

OP./C Timing Control: To display the Tj(0), Ti(z) 
functions on the TR-48 scope, both solutions are run 
in RO controlled by a * 'master timer program" on 
the DES-30. This program (Figure 5) generates a 
12 ms OP- 12 ms 10 signal by dividing the 1 KC 
square wave signal; first, a blip is produced every 
3 ms for the IC track/ store control program. On 
every fourth of these blips, a flip flop is triggered. 
The FF output is, therefore, high for 12 ms and low 
for 12 ms. 
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Figure 5: Master Timer Program 

During the first run, because both solutions are 
operating simultaneously, the track/ store units 
have not yet established the correct outputs and the 
LC.'s of the pseudo- serial solution are incorrect. 
This situation is corrected on the second run. 

Fasf-fo-S/ow Converter (FSC): Solutions generated in 
the RO mode^ can be plotted on conventional x-y 
plotters by means of a * 'fast- to- slow converter" 
program. This program consists of two track/ 
store units cascaded as shown in Figure 6a. (Tim- 
ing for this program is shown in Figure 6b.) 




T* PLOTTCIt YAKl* 



Figure 6a. Program for Fast-to-Slow Conversion (FSC) 
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Figure 6b. Timing for Fast-to-Slow Conversion 

The control signal, a, for the Fast- to- Slow Con- 
verter" could be generated with two integrators and 
an electronic comparator. To conserve integrators, 
however, as well as to demonstrate DES- 30 capabil- 
ities, the control signal is generated by an all- 
digital program. Thisprogram operates as a vari- 
able on -time mono stable timer, generating increas- 
ingly longer pulses in equal increments. These 
increments can be made so small (47 ms in this 
program) that the output of the FSC can be plotted 
directly without smoothing. 

With this control signal applied, T/S 34 **mem- 
orizes" the value of the high speed solution at times 
T, 2T, 3T,. . . nT, where T is the basic increment 
(47 ms). 

T/S 35 accepts the high speed solution at T, then 
2T, 3T, etc.; its output, therefore, is the value of 
the function at these instances, changing progres- 
sively at each RO run. Thus, after the first RO run, 
the output of T/S 35 is the value of the solution at 
t =T; after the second RO run, at t =2T; etc. If the 



tNo distinction is made between DBS' 30" control led OP-IC cycle and 
TR-48 RO-timer-control led OP'IC cycle; both situations are referred 
to as "RO" in this paper. 



increments, T, are small enough, the output of 
T/S 35 will change only slightly between each 
operate cycle. With T = 47 ms and the OP time 
being 12 ms, 256 runs are required for T/S 35 to 
**go through" the entire solution. Since 256 RO 
runs require 6.16 seconds, we obtain a function 
slow enough for plotting. 

"3-0" Ploffing: The focus of the demonstration is 
the automatic plotting of a surface" of tempera- 
ture, T, as a function of both distance, z, and time, 
6, under DES-30 logic control. The '^3-D" effects 
are achieved as follows: The function To(^) is 
plotted with a constant negative bias, -3v, on the y 
axis; Ti(6)^ with -2v on the y axis and +lv on the 
X axis; T2(0), with -Iv on the yaxis and +2v on the 
X axis; and T3(^), with no bias on y and +3v on x. 
This biasing produces the depth effect so that, for 
example, the Tj(^) functions appear to begin at equal 
intervals along a '*Z" axis (in effect an isometric 
projection). 

The functions T^(z) are manipulated as follows. 
Each is displaced by increasing and equal amounts 
along the x axis; for example, Ti(z) starts at 
+ 1.5 V with respect to To(^); T2(z) at + 3 v; 
T3(z) at + 4.5 V and T^(z) at + 6v (where To(^) 
ends). In addition, each Ti(z) is superimposed on 
a 45° bias (equal ramps on both x and y axes) to 
simulate the depth effect. Also, a constant -3v 
bias on the y axis is needed since T^{d) is plotted 
with this bias. 

The z-axis has only 3 intervals (as opposed to the 4 
intervals on the 0 axis) and because the z-axis in 
reality is a 45° line, the length of the time base sig- 
nal when plotting Tj^(z) must be shorter than the one 
used in plotting Tj(^). Since, in this example, Ti(z) 
has only moderate curvatures, this variation in the 
time base is achieved by speeding up the FSC by a 
factor of 2 and reducing the time base integration 
rate when plotting T^(z) curves. This comes about 
because the FSC is actually determining the plotting 
speed, while the integration rate of the time base 
integrator is similar in function to the plotter's 
scale factor controls. 

For the 3-D temperature surface plot, the DES-30 
must perform the following tasks: 

(1) Select, by means of a set of 8 electronic 
switches, one of the eight T plots and apply 
it to the FSC input. 

(2) Control the pen drop and lift, providing 
proper delays and mode control commands to 
assure clean plotting. 
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(3) Control the time base generation and the ap- 
plication of proper biasing for 3-D effects. 

(4) Stop the plotting after the surface'' is com- 
plete. 

A sketch of a typical **3-D" plot obtained with this 
program is shown in Figure 7. 

SCALING 
Equations 



Parallel 
Solution 
oil T = T(e) 



dT 
o 

de 

dT^ 
de" 

de 

dT. 



= 2 



T - T 
_1 o 



- 2T + T 
2 10 

(Az)^ 



(Az)' 



de 



(Azf L 



Parameiers and Constants 
Label 



(7) 



(8) 



(9) 



(10) 



Pseudo- 
Serial 
Solution 
all T = T(z) 



< 





T - T 

1 0 
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dT^ 
dz^ 






T - T 

3 2 


^ 2 
dz 


Ae 


d\^ 
^dz^ 


T - T 

4 3 

Ae 


✓ \ 





(11) 



(12) 



(13) 



(14) 




Figure,?' X-Y Plotter Display - Axes Shifting and Biasing for 
3-D Effects are Accomplished Automatically 



Meaning 



Units 



Value 



L 


1/2 slab thickness 


ft 


0o25 


c 


specific heat of slab 


Btu/lb - OF 


0.3 


P 


weight density of slab 


Ibs/ffc^ 


139 


h 


air film coefficient 


Btu/hr-ft^ - Of 


5 


k 


thermal conductivity of slab 


BtuAr-ft - OF 


0.63. 


Oi 


k/pc 


ft^Ar 


0.015 


hL 
k 




dimensionless 


2.0 






lAr 


0.24 


X 


distance along slab thickness 


ft 


0-0.25 


t 


time 


hrs 


0-6 


z 


x/L 


dimensionless 


0-1 


e 


(a/L^)t = 0.24t 


dimensionless 


0-1.44 


Az 


slab segment 


dimensionless 


1/3 


Ae 


"time" segment 


dimensionless 


0.48 


T 


computer time, non RO 


seconds 


0-6 



See Time 
Scaling Section 
on P. 8 



Unsealed Compufer Program: The parallel Solution 
is programmed as shown in Figure 9. The pseudo- 
serial solution is programmed as in Figures. The 
quantities to be scaled are thus evident. (See next 
section below) . 



o— 




dTi 



o— 





Figure 8: Typical Unsealed Block in the Pseudo Serial 
Solution 



Amplitude Scaling (Normalized) 

Parallel Solution 



Variable 



T - T 

2 1 



- T^ 

A 3 



^3 ■ ^ ^1 



Max. 

Value 



1000^ F 



lOOO^F 



lOOO^F 



1000^ F 



lOOO^F 



lOOO^F 



lOOO^F 



lOOO^F 



lOOO^F 



lOOO^F 



Computer 
Variable 



[loooj 
[loooJ 
[loooJ 

T3 

[loooj 

T^ - T^ 
L 1000 J 

[ 1000 J 

- T^ 
1-3 21 

L 1000 J 



- T^ 

L 1000 J 



T^ - 2T + T 

L 1000 J 

T - 2T + T 
I 3 2 li 



1000 



d^T3=f Az(T^-T3) 1000°F 



LioooJ 



h -1 



T.(e) 





Q tX '<1^ 



^•0 



Figure 9: The Unsealed Parallel Solution 

Pseudo -Serial Solution 



Variable 



T3(z) 



dT. 



dz 

dT^ 
2 

dz 

dT^ 
3 

dz 
dT 

4 
dz 



T - T 

1 C 



Max, Value 



lOOO^F 



1000° F 



lOOO^F 



lOOO^F 



1000 FAr 



1000 F/hr 



1000 F/hr 



1000 l/hr 



lOOO^F 



Computer Variable 
T^ 

LioooJ 

T^ 

[loooJ 

Tg 

LioooJ 

T 

1—1 

L loooj 

T^ 

L loooJ 

T^ 

LioooJ 

Tg 

LioooJ 

T^ 

LioooJ 

T - T 

L 1000 J 
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Pseudo-Serial Solution (Continued) 
Variable Max. Value Computer Variable 



T - T 

3 2 



4 3 



lOOO^F 



lOOO^F 



1000° F 



T - T 

2 li 



[ 1000 J 

L 1000 J 
L 1000 J 



The scaled equations (not repeated here) are 
identical to equations (7) through (14) except that 
a factor of 1/1000 appears on both sides of each 
equation. This is due to the 1/1000 scale factor 
which appears in both the variables and their 
derivatives. 

Time Scaling: Since Computer time represents 
time, (6), for the parallel solution integrators, 
while it must, simultaneously, represent distance, 
(z), for the pseudo- serial solution, there are two 
''time scale" factors, jSj and jgj^. The run time is 
no longer arbitrary since it represents z, and any 
change of the run\time will require rescaling. In 
addition, the constants involved in equation (1) are 
normally stated in terms of units such that time is 
measured in hours rather than seconds. These 
factors must be born in mind when trying to inter- 
pret the computer program. 

The two **time scale" factors, jSj and jSj, are de- 
termined as follows. Since r =i8j 0 where ris 
computer time in secorids, then r =0.24)5jt where t 
is real time in hours (recall 9 =0.24t). If we let 
0,24Pj =1, then 1 second computer time represents 
1 hour real time. Hence, ft =1/0.24 =4.16 seconds. 
(Note the unusual units of this time scale factor. 
This is due to the normalization of the system equa- 
tion, replacing t by Q) . 

For the pseudo-serial solution, r =iSiZ. To de- 
termine jSi, it is necessary to know how long a run 
is needed for the parallel solution and adjust jSl 
accordingly so that when z =1, ftz equals the run 
time. It is anticipated that six hours will be the 
maximum period of interest; hence, r = 6 seconds 
when z =1. Therefore, jS^ =6 seconds. (Again note 
the units of jS-) . 

Note that since there are four stations in space 
(j = 0, 1, 2, 3), Az = 1/3. In time, (9), there are 
five stations (i =0, 1, 2, 3, 4) but the computer 
doesn't know this. There are only 4 dynamic 
blocks in the pseudo- serial solution, corresponding 



to i =1, 2, 3, 4 only since for i = 0 the Ti(z) function 
is simply a constant representing the original uni- 
form temperature distribution in the slab. Hence, 
A e must be 1/3 of whatever 0 value is represented 
by 6 hours of real time, namely 0 =0.24t = (0.24) 
(6) =1.44. Hnece, A0 =(1.44) (1/3) =0.48. 



To summarize: 



jSj = 4. 16 sec 



jSj^ = 6 sec 



(1 sec, computer time 
= 1 hour real time) 



(1 sec. computer time 
= 1/6 slab thickness 
= 0.5'») 

Az =1/3 dim. less (4 stations in space) 

A0 = 0.48 dim. less (4 effective stations in 

time) 

The selection of the OP period in the RO mode 
is not arbitrary (since run time equals distance for 
the pseudo- serial solution). Since activating the 
TR-48 Time Scale bus PES-30 FAST/NORMAL 
switch in FAST) causes a 500:1 speed up in solu- 
tion rate, the OP period must be 6/500 seconds or 
12 milliseconds. 

The time scale factors are implemented, of course, 
by manipulating integrator gains rather than writing 
time-scaled equations. 



HOW THE PROGRAM WORKS 



Appendix IV shows the detailed program. It is 
divided into 3 major functional blocks: Analog 
Program and Master Timer; Y-Axis Switching and 
Fast-to-Slow Converter; and, 3-D Plot Control. 
(The entire diagram can be obtained by combining 
these three blocks pages on one sheet.) 

The system for which this program is written con- 
sists of the following: 

(1) A fully expanded EMC TR-48 (16 integra- 
tors and 32 amps.), plus two Quad ampli- 
fier groups (for a total of 36 amps.); with 
sweep and blanking inputs to the RO scope 
(allowing it to operate in all modes, not 
only RO) ; pen drop input (logic 1 = + 5v = 
drop pen); four separate y-axis to the RO 
scope and x, y inputs to the plotter (y^ scope 
input not shared with plotter y input) ; and 
with 8 electronic comparator units (8 com- 
parator networks, 16 DA switches, 16 track/ 
store capacitor networks) . 
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(2) A fully expanded DES-30 with 40 flip flops, 
48 AND gates (2-input), 4 down counters, 4 
MT-DIF groups, 18 DA trunks, 7 AD trunks. 

(3) 1110 plotter 

(4) External wide-band oscilloscope for setting 
the MT and checking the logic program dy- 
namically. 

Analog Program: The parallel solution consists of 
integrators 36, 37, 42, 43; amplifiers 24, 25, 40, 
41, 28, 29, 32, 33; and pots 45 through 49 and 59. 

Pot 04 (with amplifiers 16 and 17) establishes the 
initial uniform slab temperature Tq(z) as identical 
IC's for the integrators of the parallel solution. 
Thus a convenient control for this parameter is 
provided. 

Pot 49 contains both h (air film coefficient) and k 
(thermal conductivity of slab) and thus provides 
control over both. 

Pot 59 has on it the cooling air temperature T^. 

The pseudo- serial solution consists of integrators 
00, 01, 12, 13, 30, 31, 44, 71; amplifiers 04, 02, 15, 
19, 47; and pots 03, 00, 16, 15, 23, 52, 50, 53. Note 
that P04 provides the fimction Tq(z) (a constant) 
which is needed in the first block of the pseudo- ser- 
ial solution. 

The initial conditions for the Ti(z) integrators are 
zero because of the symmetry considerations al- 
ready explained. The initial conditions for the 
Tj(z) integrators are obtained by the parallel-to- 
serial IC pickup program which consists of track/ 
store units, 26, 27, 38, and 65; and GPR 26, AND 
gates 16A through 16E, and MT6. The function 
of this program is to sample (and hold) the value 
of the temperature as a function of 9 (time) at four 
equally- spaced instances of the OP period (cor- 
responding to i =1, 2, 3, 4 but not i =0). 

These values are the correct I.C.'s for integra- 
tors 01, 13, 31, 71 respectively. This is achieved 
as follows: The OP period, which is 12 milli- 
seconds long, is indicated by flip flop 270 being 
high (this FF is part of the master timer which 
will be explained shortly) . When this is the case, 
the output of AND gate 16E is four equally- spaced 
(3 ms apart) blips; the first blip occurs 3 ms after 
going into OP. These blips are derived from the 
master timer program. Each time such a blip oc- 
curs, MT6 generates a 100 microsecond pulse. 
This pulse is applied to the TRACK (**DIG") input 
of the track/ store units in sequence, beginningwith 



TS-26 down to TS-65. The sequencing is controlled 
by GPR-26, which is patched as a ring shift regis- 
ter. During the IC period (FF 270 low), a bit is 
forced into stage D of this GPR through the L-E 
control. Thereafter, each blip that sets MT6also 
shifts this bit one stage ahead. The present posi- 
tion of this bit determines which of the four AND 
gates 16A through 16D will transmit the 100 micro- 
second signal generated by MT6. 

r— — a \ I I 

MT -H- H- n n n n 
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Figure 10: Timing of the Parallel-to-Serial IC Pickup Control 

Amplifier 66 is used as a buffer. The low input 
impedance of the track/ store units with the small 
capacitor causes ringing if the driving amplifier 
drives 3 or more track/ store units simultaneously. 
Note that amplifiers 36 and 66 drive only two T/S 
units each. 

Moster Timer: The function of the master timer is 
to generate a 12 ms OP - 12 ms IC signals as well 
as the four equally spaced blips during OP for the 
parallel-to- serial IC pickup just described. This 
master timer program consists of down counters 
3 and 2, flip flops 27B and 27C, and AND gates 16F 
and 17F, The function of MT-7 and DIF-7 will be 
explained shortly; for the time being assume they 
are removed (Cq of DC-2 to Trigger of FF-27C) . 

AND gate 16-F and FF-27B buffer (clock) the asyn- 
chronous IKC square wave source. The output of 
DIF-6 is a blip every 1 millisecond. DC-3 preset 
switch is set to 02, so its Cq output is a blip every 
3 ms (recall that one carry in is required to load 
the DC with the preset number, during which time 
EC^ is disabled, so the total number of carry-ins 
per carry-out is the preset number plus one ) . The 
Cq of DC-3 is thus the signal needed for the 
parallel-to- serial IC pickup control. 

DC-2 counts 3 + 1 =4 of these 3 ms blips, so that 
its Cq is a blip every 12 ms. This blip may be used 
to trigger FF-27C so that it is high for 12 ms and 
low for 12 ms. The outputs of FF-27C can then 
drive the OP-RST inputs which control the mode 
of the TR-48. 

Two modifications were found necessary to the last 
paragraph: first, the beginning of the OP signal 
had to be delayed from the Co of DC-2 by some 
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60 microsecond. The reason for this delay will 
be explained under the fast- to- slow converter 
description. The delay is achieved by setting MT- 7 
from the of DC-2 and trailing edge differentia- 
ting (DIF-7) the output of MT-7. Thus the output 
of DIF-7 is a blip occurring about 60 jUS after the 
Cq of DC-2, and it's this differentiator output that 
is used to trigger FF-27C. Second, in order to 
assure that the last TRACK signal (the one for TS- 
65) in the IC pickup program does not extend beyond 
the OP period, the output of MT-6 is OR'ed in OR 
gate 18F with the output of FF-27C, and the output 
of this OR gate is the OP-IC signal. The effect 
of this modification is felt only at the end of the OP 
period; rather than going to IC as soon as FF-27C 
resets, the TR-48 waits for MT-6 to complete its 
100 Ids pulse, so that track/store 65 tracks only 
within the OP period. This is best explained by 
the timing diagram below. (Figure 11). 



-^c. I \ \ \ L I I 

DC Z C, I I 

| l J 

^•'c ;| I 

i | -n*-*"^*n n fji 

Figure 11: Master Timer Timing Diagram 

AND gates 18D and 18E allow dropping both OP 
and RST inputs to logical zero for the HOLD mode. 
The TR-48 is thrown into HOLD whenever the plot- 
ter pen is about to be lifted, to prevent streaks. 
The pen lift operation will be discussed in the 3- 
D Plot Control section. 

The signal out of AND gate 17F is a blip occurring 
about 60 microseconds prior to going into OP. 
This signal is labeled MTP (Master Timing Pulse) 
and is used in the fast- to- slow converter . 

Y Axis Switching: Since the plotter will be required 
to plot 8 different functions (four parallel and four 
pseudo- serial solutions), these 8 plots must be 
switched into the y axis of the plotter in some 
sequence. (Fast- to- slow conversion must, of 
course, take place first; this will be described 
shortly). DA switches 14, 15, 22, 23, 02, 03, 46, 47 
accomplish this y- axi s- switching . Each switch con- 
trols one of the 8 functions and conducts when the 
signal on the DA trunk patched into its DIG (con- 
trol) input is high (= +5v). The signals on DA 
trunks 0-5, 18, 19 are derived from the 3-D Plot 
Control which will be explained later. 

The switches are patched in a somewhat un- 
orthodox manner (which may well be the standard 
configuration in TR-48/DES-30 systems). The 



input to the switch with its lOK bult-in input re- 
sistor is not being used at all, since with this low 
resistor value the switch is non-linear and must be 
placed in the feedback loop of an amplifier, as 
shown in the Electronic Comparator Manual; but this 
implies SPDT operation requiring 2 switches per 
signal to prevent a no-feedback condition when the 
switch is non-conducting. Instead, the GJ (Gate 
Junction) input is used with external lOOK resistors. 
This, incidentally, permits summing into the switch, 
as shown on switches 14, 15 and 16. The constant 
voltages summed into these switches in addition to 
the functions Tq (^), T^ (6), T2 (6), are bias 
voltages for 3-D effects. These will be explained 
in the 3-D Plot Control section. The outputs of the 
switches ('*B'') are all patched into the SJ of A23, 
whose output is minus whichever signal is per- 
mitted to pass through any conducting DA switch. 
The '*0" point of the DA switches (patched to the 
amplifier output in SPDT configuration) is not used. 

FasMo-Slow Converter: The function of the fast-to- 
slow converter (FSC) is to convert solution wave- 
shapes produced in RO speeds to identical wave- 
shapes produced at rates slow enough for plotting. 
Track/ store units 24 and 25 perform this conver- 
sion. A control signal (the true output of which is 
patched to TRACK of TS-34 and its complement to 
TRACK of TS-35) is generated by a '*variable-ON- 
time monostable'* (VOTM) program consisting of 
GPR 21, 14, 20, 15, down counters 0 and 1, and 
associated flip flop and gates. 

The principle of operation of the FSC is this: if 
the control signal is high from the beginning of the 
OP period until time tQ (where to < OP period, 
12 ms in this program), then T/S 34 stores the 
value of the solution at tQ. T/S 35 begins tracking 
only when T/S 34 is storing, so its output is also 
the value of the solution at tg. As long as the 
control signal remains unchanged, the output of 
T/S 35 remains a constant (this is not true of T/S 
34). If the time t^ is slowly increased (the control 
signal remains high for longer portions of the OP 
period), then the output of T/S 35 will slowly trace 
the values of the high speed solution. 

The three major components of the control signal 
generation program (VOTM) are: an 8-bit binary 
down counter, consisting of GPR 21, 14 and gates 
12F, 12E, 13E, 13F, llA, and lOA; an 8-bit binary 
up counter, consisting of GPR 15 and GPR 20; and 
a two-decade BCD down counter pC 0) . Note that 
GPR 15/20 is also patched as a ring shift register, 
and that GPR 21/14 can accept serial inputs. 

DC 0 and FF-29D operate as follows. Every time 
FF-29D is set, DC 0 cotmts down for 47 ns; then 
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its carry-out decrements the binary down counter 
GPR 21/14 and also tries to reset FF-29D through 
AND gate 12C. As long as GPR 21/14 contain a 
number other than zero, FF 29D will not reset, 
so DC 0 reloads itself with 46 and starts counting 
down again. After GPR 2 1/14 has been decremented 
to zero, the Cq of DC 0 resets FF29D and stops 
the counting. The output of FF 29D thus remains 
high for a period of 47 (N + 1) microseconds, 
where N is the number stored in GPR 21/14 when 
FF 29D set set. 

This number N is loaded into GPR 21/14 serially, 
before FF 29D is set, from the binary up counter 
GPR 15/20. This is accomplished as follows. 
Whenever FF 27D is set, DC 1 counts down 8 
clocks and then resets FF 27D. Thus FF 27D out- 
put is high for exactly 8 clock periods, and this 
signal is used to shift the contents of GPR 15/20 
into GPR 21/14. The Cq of DC 1 (which resets 
FF 27D and thus stops the shifting) also sets FF 
29D which starts the process described in the last 
paragraph. 

The signal that sets FF 27D (thus starting the 8-bit 
shifting) is either the output of FF 28A or the 
output of AND gate 19B or IOC. Assume, for the 
time being, that this signal comes from AND gate 
IOC, in which case it is the MTP blip from the 
master timer . (Recall that the MTP indicates the 
beginning of the OP period, and that it occurs about 
60 microseconds prior to the start of the OP 
period) . 

When this is the case, then the sequence of events 
that take place is as follows: 

(1) MTP 

(2) GPR 15/20 shift 8 zeros into GPR 21/14 

(3) 47 MS microsecond pulse at output of 
FF 29D 

(4) Wait for next MTP. 

Thus, the output of the FSC (T/S 35) is the value 
of the high speed solution at about 10 microseconds 
before going into OP; namely this is the IC value. 

If, before the next MTP, GPR 15/20 is incremented 
(counts up) once, then the sequence above will 
repeat, except step 2 will read: shift 7 zeros and 
a 1; and step 3 will read: 94 microsecond pulse, 
etc. If GPR 15/20 counts up once before every 
MTP, then step 2 in the above sequence will read: 
shift an 8 bit binary number N into GPR 21/14; 
and step 3 will read: 47 (N + 1) microsecond pulse 
etc. 



The incrementing of GPR 15/20 is controlled by 
FF 27A, 28A and gates 17E, 18E, andl9E. Assume 
for the moment that the MTP is routed to AND gate 
lOB (rather than IOC) . The first MTP starts FF 
27D - DC 1 as before, this time through AND gate 
19B, but it also sets FF 27A. When the next MTP 
arrives, it does not pass through AG 19B because 
FF 27D is set; instead it goes on through AG 17E 
to set FF 28 A and also to increment GPR 15/20 
once through OR gate 19E. 

If the signal is high, which indicates that the solu- 
tion being plotted is one of the pseudo- serial 
functions, then the output of OG 19E will be high 
for 2 clock periods, so it will increment GPR 15/20 
twice. In this way the FSC is made to increment 
twice as fast as before; this requirement satisfies 
(partially) the need to snorten the time base of the 
pseudo- serial solutions (see 3-D Plot Control 
explanation) . Note that whether GPR 15/20 is incre- 
mented once or twice, the signal which starts 
FF 27D and the 8-bit shift (dumping GPR 15/20 
contents into GPR 21/14) is the output of FF 28A, 
so the shifting will not start until the counting 
(incrementing) is completed. 

The signal that controls the routing of the MTP to 
AG IOC (* 'generate FSC control, but do not incre- 
ment GPR 15/20") is the PEN DROP signal. This 
signal, which indicates that the pen is down and 
ready to plot, is generated by the 3-D Plot Control 
section. As long as the pen is up, the MTP is 
routed to AG IOC, and the output of T/S 35 is the 
IC value. When the pen is down, the MTP goes 
through AG lOB, GPR 15/20 is incremented on 
every MTP, and the output of T/S 35 slowly traces 
the high speed solution. 

When GPR 15/20 is full (reading 255), the carry- 
in from OG 19E causes a carry-out; this Cq is the 
FSC FINish signal, indicating that 1 complete plot 
has been traced by T/S 35. 

The HOLD signal OR'ed into the FSC control sig- 
nal (OG 17A) is generated by the 3-D Plot Control. 
It indicates that the pen is in the process of lifting 
at the end of the plot; T/S 35 is therefore asked to 
store (HOLD) so that there are no ** streaks" on 
the paper. 

3-D Plot Control: The major components of the 3-D 
Plot Control program are: y- axis switching control 
(GPR/22/23, comparator 11 and 03); logic-con- 
trolled time base (integrator 20, DA switch 11, 
pots 10 and 20, OR gates IIC/IID and IIB; x-axis 
bias (amplifier 22 and the associated pots); pen 
drop control and miscellaneous items (rest of 
diagram). 
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GPR 22/23 are patched as a ring shift register. 
Comparator 03 supervises" the initialization of 
GPR 22/23 by loading a bit in stage A on the first 
DES clock, after which the comparator '*with- 
draws". The compairator output is low in PS, so 
that the complementary outputof ADtrunkS is high, 
and since the DES is put in RUN before the TR-48 
is put in SL, the loading of stage A occurs. Then 
the TR-48 is put in SL, its mode is cycled between 
OP and IC by the master timer , comparator 03 out- 
put is high and AD trunk 8 complementary output 
is low. The overall effect is as if this situation 
obtained: 



The outputs of GPR 22/23 are the signals which 
control the DA switches in the Y-Axis Switching 
program previously described. Each time a plot of 
one function is finished, the EOP (End of One Plot) 
signal (which will be explained shortly) steps GPR 
22/23 and thus selects the next function to be 
plotted at the input to the FSC (output of A23). 

The outputs of GPR 23 are OR'ed in OG IIE/IIF 
to generate the signal 8 which signifies that the 
pseudo- serial solutions are being plotted. This 
signal is used in the VOTMprogram. It also causes 
DA switch 10 to conduct. This applies the output of 
integrator 20 to the y-axis to achieve the 45°-line 
(z-axis) effect, since when switch 10 is conducting, 
the output of integrator 20 is applied simultaneously 
to the X and y axes. 

Comparator 11 simply acts as an inverter saving 
a DA trunk which would otherwise be needed to 
send 8 to the TR-48. When 5 =1, the sum of the 
comparator inputs is +4 volts (pot 14 is set to 
0.1) so its complementary output is low; when 
5=0, the comparator input is -1 volt and its com- 
plementary output is high. Thus DA switch 11 is 
non-conducting when 8 =1, reducing the time base 
integration rate when the pseudo- serial solutions 
are plotted as explained previously. 

The outputs of GPR 22/23 are also used to generate 
the X-axis bias for each plot. However, since the 
outputs of the DA trunks are only approximately 
+ 5v when the input is logical 1 and only approxi- 
mately 0 volts when the input is zero, the pots (25, 
26, 27, 28, 40, 43, 42) must be set so that the out- 
put of A22 gives the correct bias. In addition P29 



is added so that when all switches are off, it can 
be adjusted to make A22 output actually zero volts. 

The plotting process is started by pushing MP 4C. 
This sets FF 29C, whose complementary output, 
the STOP signal, becomes low. FF 28B is also set 
at the same time, causing the plotter pen to drop. 
Integrator 20 starts integrating (provide x-axis 
time base), since the output of OG IIB is low and 
the complementary output of OG IIC/IID is high 
(OP input = 1, IC input = 0 on integrator 20) . 

The plotting proceeds until the FSC FIN signal 
arrives, indicating the completion of one plot. This 
resets FF 28B (pen up) and sets FF 28C (HOLD 
mode) . 

The HOLD mode is retained for about 1 second as 
follows. FF 29A is enabled by the differentiated 1 
CPS square wave. When FF 28C goes high, FF 29A 
waits for the next 1 CPS blip from DIF 9, and then 
it goes high until the next 1 CPS blip. FF 28C 
resets on the same blip that resets FF 29A, so 
that its output is guaranteed to be high for at least 
1 second, and possibly up to 2 seconds (this latter 
situation occurs when the output of 28C goes high 
just after 1 CPS blip). 

The blip that resets the HOLD flip flop also goes 
through AG 19A as the EOP signal — End of One 
Plot — which steps GPR 22/23 and also resets the 
VOTM program so it's ready for the next plot. In 
addition EOP goes through AG 18C to set MT 9 
which generates a 0.5 second (appx.) signal. This 
signal throws A20 into IC. At the end of the 0.5 
second delay, trailing edge differentiator 8 sets 
FF 28B (pen down) again, and the process repeats. 

If LP 5B is depressed, the complementary output 
of AG 17B is low. This prevents FF 29C from being 
reset, so that the plotting process will continue 
indefinitely. If LP 5B is released, then when the 
bit in GPR 22/23 returns to stage A after completing 
8 plots, it causes the next EOP signal to reset 
FF 29C rather than set MT 9, thus stopping the 
process. 
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APPENDIX I: BOUNDARY EQUATIONS 



The ''parallel" equation at the center of the slab is dTg , ^ 

dTp ^ Ti-2To+T.i 

But T^i =Ti by symmetry, hence by replacing -r — by -— ; then 

o z ^ZX ^ 



The equation at the air-glass interface is and 

dTg ^ T4-2T3 + T2 dTg 

de (Az)^ 

(10) 

where T4 is an imaginary point outside the slab. 
T4 can be obtained from the heat transfer equation 
across the glass-air interface: 
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APPENDIX II: STATIC TEST 




"Static Test" Copy 
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